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Abstract. Calculations of the one-hole spectral function of 16 O for small missing 
energies are reviewed. The self-consistent Green's function approach is employed 
together with the Faddeev equations technique in order to study the coupling of both 
particle-particle and particle-hole phonons to the single-particle motion. The results 
indicate that the characteristics of hole fragmentation are related to the low-lying states 
of 16 and an improvement of the description of this spectrum, beyond the random 
phase approximation, is required to understand the experimental strength distribution. 
A first calculation in this direction that accounts for two-phonon states is discussed. 
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1. Introduction 

Recent advances in nuclear theory have generated accurate predictions for the spectrum 
of most p shell nuclei (see, for instance, Refs. PUI2])- ^ ^he same time, other techniques 
are becoming available to describe larger systems and to account for the effects of 
the continuum [HIE]- For medium and heavy systems, relevant information regarding 
correlations has been obtained by studying the nuclear spectral function. This was 
mainly done by means of variational calculations [H] and the self-consistent Green's 
function (SCGF) |jj approaches. 

Once a given nuclear Hamiltonian is chosen, quantities such as the spectroscopic 
factors are defined uniquely in term of the exact solutions of the many-body problem. 
Thus, their knowledge gives direct information on the correlations induced by that 
specific nuclear force. Experimentally, (e, e'p) reactions have provided results for knock 
out from orbits both close to [7J |S| and far from [5] the Fermi energy. Although a 
consistent calculation (based on the same Hamiltonian) of the initial and final states has 
so far been possible only for specific cases jTUj, several analyses [TTJ[T2J[TH1[I| suggest that 
the experimental cross section can be described by standard phenomenological realistic 
interactions. This leads to fragments at small missing energies that have spectroscopic 
factors of about 60-70% [7j. Moreover, recent measurements of the spectral function 
at high missing energies and momenta ^3] appear to be consistent with the tail due 
to the short-range and tensor correlations (SRC) that are induced by nuclear forces 
having a repulsive core. Whether (and how) softer NN interactions can describe these 
measurements is an open (and interesting) question. 

For the case of 16 O, there still exists a substantial disagreement between the 
quenching of spectroscopic factors extracted from the experiment (HI CHI E2] and 
theory [TI3 HH1 El- The latter results suggest that the reasons for this discrepancy 
should be looked for in the effects of long-range correlations (LRC) and in particular in 
the couplings of single-particle (sp) motion to low-energy collective excitations. In this 
contribution we report about the work done along this line in Refs. [TTf IT8] in order 
to tackle the above issues for 16 0. Sec. El describes the SCGF and Faddeev formalism 
employed to couple sp and collective phonons [15] . The results for the hole spectral 
function are discussed in Sec. El These calculations show that a proper description 
of the experimental spectral strength requires an improvement of the spectrum which 
goes beyond the random phase approximation (RPA). A first step in this direction that 
includes the propagation two-phonon states is reported in Sec. HI 

2. Faddeev approach for the single-particle Green's function 

We consider the calculation of the sp Green's function 
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Figure 1. Diagrams contributing to the irreducible self-energy £*. The double lines 
represent a dressed propagator and the wavy lines correspond to a G-matrix (that is 
used in this work as an effective interaction). The first term is the Brueckner-Hartree- 
Fock potential while the others represent the 2plh/2hlp or higher contributions that 
are approximated through the Faddeev TDA/RPA equations. 



from which both the one-hole and one-particle spectral functions, for the removal and 
addition of a nucleon, can be extracted. In Eq. (JTJ, X£ = (^^ +1 |c^,|^o) O^a = 
( 1 I r ^ _1 |c a | 1 i r Q )) are the spectroscopic amplitudes for the excited states of a system with 
A + 1 (A - 1) particles and the poles e+ = E^ +1 - E$ = E$ - E^~ l ) correspond to 
the excitation energies with respect to the A-body ground state. The one-body Green's 
function can be computed by solving the Dyson equation 

9afi(w) = g°ap(u) + ^g%(u)YZtf(u)g S {}(u) , (2) 

where the irreducible self-energy T,* s (u) acts as an effective, energy-dependent, 
potential. The latter can be expanded in a Feynman-Dyson series |2lH |2*T] in terms 
the exact propagator g a p{u), which itself is a solution of Eq. (J2J). In this expansion, 
T,* 5 (u) can be represented as shown in Fig. [U by the sum of a dressed Hartree-Fock 
potential and terms that describe the coupling between the sp motion and more complex 
excitations i>j. It is at the level of the 2plh/2hlp propagator, R(u), that the correlations 
involving interactions between different collective modes have to be included. 

The SCGF approach can be initiated by solving the self-energy and the Dyson 
Eq. (j2J) in terms of an unperturbed propagator g%(u>). The (dressed) solution g a p{oj) 
is then used to evaluate an improved self-energy, which then contains the effects of 
fragmentation. The whole procedure is iterated until self-consistency is reached. Baym 
and Kadanoff showed that a self-consistent solution of the above equation guarantees 
the fulfillment of the principal conservation laws 



2.1. Faddeev approach to the self-energy 

In the following we are interested in describing the coupling of sp motion to ph and 
pp(hh) collective excitations of the system. All the relevant information regarding the 
latters are included in the Lehmann representations of the polarization propagator 

(«tfl4«*l*»> <*£l4*l*o> 



n a/3 , 7 5(w) = 
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Figure 2. Example of a diagram appearing in the 
all-orders summation generated by the set of Faddeev 
equations. 
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which describe the excited states of the systems with A and A±2 particles, respectively. 
In general, Eqs. (jUJ) and (j3J) are the exact solutions of their respective Bethe-Salpeter 
equations (BSE). In the calculation of Sec. El these have been approximated by solving 
the dressed Tamm-Dancoff/RPA (DTDA/DRPA) equations (2H1 Ell, which account for 
the effects of the strength distribution of the particle and hole fragments. The inclusion 
of correlations beyond RPA is considered in Sec. 0] 

The ph (jSJ) and pp(hh) (@J) propagators are inserted in the nuclear self-energy by 
solving a set of Faddeev equations [2H] for the 2plh and 2hlp propagators of Fig.HJ The 
details of this approach are given in Ref. jTHj. For the present discussion it is sufficient 
to note that the motion of three-quasiparticle excitations is approached in the same way 
it is normally done for the three-body problem. Collective excitations are coupled to sp 
propagators generating an infinite series of diagrams, including the one shown in Fig. |2l 
This allows to account completely for Pauli correlations at the 2plh/2hlp level. 

3. Results for the single-particle spectral function of 16 

In the calculations described below, the Dyson equation was solved in a model space 
consisting of harmonic oscillator sp states. An oscillator parameter b = 1.76 fm was 
chosen (corresponding to hu = 13.4 MeV) and all the first four major shells (from Is 
to 2pl/) plus the l<7g/2 where included. The results of Refs. j2HIHl, suggest that this 
model space is large enough to properly account for the low-energy collective states 
if fragmentaion is accounted for. Inside the model space, a Brueckner G-matrix 
derived from the Bonn-C potential ,27] was used as an effective interaction. The 
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Shell 


TDA 


RPA 


1st itr. 


2nd itr. 


3rd itr. 


4th itr. 


Z Pl/2 


0.775 


0.745 


0.775 


0.777 


0.774 


0.776 


Z P3/2 


0.766 


0.725 


0.725 


0.727 


0.722 


0.724 






0.015 


0.027 


0.026 


0.026 



Table 1. Hole spectroscopic factors (Z a ) for knockout of a £ = 1 proton from 16 0. 
The columns 'TDA' and 'RPA' refer to the initial undressed calculations, while the 
remaining colums resulted from the first four iterations of the DRPA equations. Note 
that the iterated resulted were obtained by constraining the lowest + solution for 16 O 
at its experimantal value, which is at the origin of the fragmentation of the p 3 / 2 peak. 



short-range core of this NN interaction induce an additional 10 % reduction in the 
spectroscopic factors [Hj, which is accounted for in the solution of the Dyson equation 
through the energy dependence of the G-matrix. 

3.1. Effects of RPA correlations and fragmentation 

For an unperturbed initial propagator the TDA calculation is equivalent to the one of 
Ref. [15] and yields spectroscopic factors equal to 0.775 and 0.766 for the main py 2 
and p 3 /2 qhasihole peaks, respectively. These results are reported in Table [T] The 
introduction of RPA correlations reduces these values and brings them down to 0.745 
and 0.725, respectively. This shows that collectivity beyond the TDA level is relevant 
to explain the quenching of spectroscopic factors. We note that due to center-of-mass 
effects, the above quantities might need to be increased by about 7% before they are 
compared with the experiment j28j . 

The RPA results were then iterated a few times to study the effects of fragmentation. 
Since only the low-energy excitations are of interest here, it is sufficient to keep track 
only of the largets fragments that appear — close to the Fermi energy — in the (dressed) 
sp propagator, Eq. ((l}, while the residual strength is collected in a effective pole |29| [TTj . 
Only a few iterations were required to reach convergence. The effect of including 
fragmentation in the construction of the RPA phonons is to increase the strength of 
the main hole peaks. The pi/ 2 strength increases from the 0.745, obtained with the 
undressed input, to 0.776. Analogously the total strength in the p 3 / 2 peak rises to 
0.750. This behavior is due to the competing effect of the redistribution of the strength, 
which tends to screen the nuclear interaction. 

The converged distribution of one- hole strength is shown in the mid panels of Fig. El 
where it is compared to the experiment (top panels). The latter is characterized by 
additional small fragments close to the Fermi energy, some of which will be discussed in 
Sec. 13.21 We note that similar results are obtained for the particle strength, including 
large peaks near the Fermi level and a fragmented distribution at larger energies. This 
self-energy at positive energies has been employed recently in studying low energy 
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Figure 3. One-proton removal strength as a function of the hole sp energy eV = 
— for 16 O for angular momentum I — 1 (left) and I = 0, 2 (right). For the 

positive parity states, the solid bars correspond to results for d 5 / 2 and c? 3 / 2 orbitals, 
while the thick lines refer to Si/2- The top panels show the experimental values taken 
from [H]. The mid panels give the theoretical results for the self-consitent spectral 
function. The bottom panels show the results obtained by repeating the 3 rd iteration 
with a modified ph-DRPA spectrum, in which the lowest eigenstates have been shifted 
to the corresponding experimental values. 



proton- nucleus scattering |3U] . 

3.2. Role of the lowest excited states in 16 O 

A deeper insight into the mechanisms that generate the fragmentation pattern can be 
gained by investigating directly the connection between the spectral function and some 
specific collective states. To clarify this point we repeated the above calculations of the 
sp propagator by shifting, at each iteration, the solution for the lowest + excitation 
in 16 O to its experimental energy. The difference with respect to the preceding results 
is the appearance of a second smaller p 3 / 2 fragment at -26.3 MeV, which might be 
interpreted as one of the fragments seen experimentally at slightly higher energy This 
solution arises in the first two iterations and converges to a spectroscopic factor of 2.6%, 
as seen in Table [0 The associated p hole spectral function is shown in the lower-left 
panel of Fig. 01 This result can be interpreted by considering the fragments as 
generated by holes in the ground state and an excited + level of the 16 O core. If the 
two levels are close enough in energy, the two configurations mix together with the result 
of fragmenting the strength over more than one peak. 

The other two low-lying states of 16 O that may be of some relevance are the isoscalar 
1~ and 3~, which are reproduced by RPA type calculations at ~3 MeV above the 
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Figure 4. Examples of con- 
tributions involving the coupling 
of two independent ph phonons. 
All the diagrams of this type, ob- 
tained by considering all the pos- 
sible couplings to a ph state, are 
included in the BSE kernel by the 
two-phonon ERPA equations of 
Ref. UH. 



experiment (see Fig. El below). The lower- right panel of Fig. El shows the results for the 
even parity spectral functions that are obtained when both the 3~ and the 1" ph-DRPA 
solutions are shifted to match their experimental values. In this d 5 /2 hole peak 

is obtained at a missing energy of -17.7 MeV, in agreement with experiment. 

4. Two-photon contributions to the spectrum of 16 

The above results suggest that an improvement of the solution for the spectral strength 
would require a better description of the excitation spectrum. One important problem of 
(D)RPA is the appearance of at most one collective phonon for a given J n , T combination 
while several low-lying isoscalar + and 2 + excited states are observed at low energy in 
16 0, as well as additional 3~ and 1~ states. A good description of the spectrum of 16 
was obtained in Ref. [HI] by coupling up to four different phonons with negative parity 
(3~ and l - ). However, the self-consistent role of coupling positive parity states and 
the dressing of sp propagators were not investigated. These effects allow for the partial 
inclusion of configurations beyond 2p2h already at the two-phonon level. Moreover, the 
inclusion of two-phonon excitations represent the first correction to the DRPA equation 
generated by the Baym-Kadanoff formalism [32] • These conisist of diagrams like the ones 
of Fig. [31 that have been included in to the kernel of the BSE. The relative formalism 
has been presented in Ref. [T%] . where it is referred to as "two-phonon extended RPA 
(ERPA)". In this work the ph-DRPA equation has been solved first, using the self- 
consistent sp propagator derived in Sec. El The lowest DRPA solutions for both the + , 
3~ and channels were shifted down to their relative experimental energies and then 
they were employed to generate the two-phonon contributions for the ERPA calculation. 

We note that the solution for the first isoscalar + state in DRPA is found much 
higher in energy at ~17 MeV and it has a sharp ph character. Therefore it cannot be 
identified with the experimental state, whose shell model structure is dominated by 
4hu configurations |HH]- On the other hand inelastic electron scattering experiments 
clearly excite this state [H3|. The one-body response is described by the polarization 
propagator, Eq. Q, and therefore the total experimental strength must be represented 
by Z ni! , Eq. (J5J). This indicates a strong coupling to ph configurations (where "ph" 
actually means "quasiparticle-quasihole" , with bare np-nh configurations implicitly 
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Figure 5. Results for the DRPA 
and the two-phonon ERPA spec- 
tra of 16 O obtained using the 
dressed input propagator com- 
puted in Sec. El middle and right 
plots, respectively. In solving 
the ERPA equation, the lowest 
3 _ , l - and + levels of the 
DRPA propagator where shifted 
to their experimental energies. 
All other DRPA solutions were 
left unchanged. The ERPA so- 
lutions indicated by dashed lines 
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phonon chatacter. The experi- 
mental spectrum is shown on the 
left. 



included by the dressing of the sp propagator). On the basis of this similarity — and 
as long as one keeps in mind its limitations — it still appears interesting to shift the 
lowest + ph-DRPA solution down in energy and investigate how it mixes with other 
configurations. 

The ERPA spectrum obtained for 16 is displayed in Fig. Eland Table El together 
with the total ph strength Z nir of each state, 

z n , = £ \«\4»cp\H)\ 2 , (5) 

a/3 

and the relative occupations of ph and two-phonon admixtures in its wave function. 
An isoscalar + state with a predominant ph character is still found at ~17 MeV, 
as in DRPA, but it is now characterized by a partial contribution from two-phonon 
configurations. Table El shows that this mixing results in a lower solution at ~11 MeV, 
which is predominately a two-phonon state. In both cases the relevant configuration 
comes from the coupling of two + phonons themselves. Of course higher configurations, 
including three- and four-phonon states, should be included to reach a complete 
understanding of the + spectrum. A study of these will be pursued in the future. 

The low-lying 3~ and 1~ states are only slightly affected by two-phonon 
contributions and remain substantially above the experimental energy at 9.23 and 
10.90 MeV. However, the coupling of these to the 0^ level reproduce the second excited 
states for the same angular momentum and parity. The two-phonon ERPA approach 
also generates a triplet of states at about 12 MeV with quantum numbers + , 2 + and 
4 + . The solutions for this triplet are almost exclusively made of 3~ (g> 3~ configurations 
and therefore have a 2p2h character. A similar triplet is found experimentally at 12.05, 
11.52 and 11.10 MeV, which correspond to twice the experimental energy of the first 3~ 
phonon. 
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Table 2. Excitation energy and total spectral strengths obtained for the principal 
solutions of DRPA and two-phonon ERPA equations, including the total contributions 
of ph and two-phonon configurations to the ERPA solutions. The individual 
contributions of the relevant two-phonon states are also indicated. 



5. Conclusions 

Self-consistent Green's function theory has been applied to study 2hlp correlations at 
small missing energies for the nucleus of 16 O. The method of the Faddeev equations 
allows to treat the coupling of ph and pp(hh) collective modes to the sp motion, The 
effects of fragmentation have been included through the dressing of the sp propagator. 
This approach allows to identify the important role played by the low-lying excited states 
of 16 0. These are essential to generate many of the fragments with small spectroscopic 
factors that are seen experimentally, examples of which are the d^/2 and p%/2 states of 
15 N at 5.20 MeV and ~9 MeV. 

The main impediment in obtaining a good theoretical description of the single 
particle spectral function of 16 has been identified in the poor description of the 
excitation spectrum, as obtained by solving the standard (D)RPA equations. We 
have improved on this by computing the effects of mixing of ph states with two- 
phonon configurations. The results show that these contributions explain the formation 
of several excited states observed at low energy which are not obtained by RPA 
calculations. However, it appears that a full solution of the spectrum of 16 O with 
this method requires to consider up to four-phonon states and the interaction in the pp 
and hh channels p?T] . 
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